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Abstract 

Theoretical results for two neutrinos in the nuclei 96 Ru, 102 Pd, 106 Cd and 108 Cd 
are presented. The study employs the Hartree-Fock-Bogoliubov model to obtain the 
wave functions of the parent and daughter nuclei, in conjunction with the summation 
method to estimate the double beta decay nuclear matrix elements. The reliability 
of the intrinsic wave functions of 96 ' 102 Ru, 96 Mo, 102,i06,i08 pd and i06,i08 Cd nuclei 
are tested by comparing the theoretically calculated spectroscopic properties with 
the available experimental data. Calculated half-lives T 2 % of 96 Ru, 102 Pd, 106 Cd 
and 108 Cd nuclei for 2v (3+ (3 + , 2v (3+EC and 2u ECEC modes are presented. The 
effect of deformation on the nuclear transition matrix element M 2v is also studied. 
PACS Numbers: 23.40.Hc, 21.60.Jz, 23.20.-g, 27.60.+j 



1 Introduction 

The two neutrino double beta (2u (3(3) decay and the neutrinoless double beta (GV (3(3) 
decay can occur in four different processes: double electron {(3~ (3~) emission, double 
positron {(3 + (3 + ) emission, electron-positron conversion {(3 + EC) and double electron cap- 
ture {ECEC). The later three processes are energetically competing and we shall refer to 
them as positron double beta decay (e + DBD) modes. The 2v (3~ (3~ decay is allowed in the 
standard model of electroweak unification (SM) and the half-life of this process has been 
already measured for about ten nuclei out of 35 possible candidates. Hence, the absolute 
values of the nuclear transition matrix elements (NTMEs) M^ v can be extracted directly. 
Consequently, the validity of different models employed for nuclear structure calculations 
can be tested by calculating the M 2u - In case of 2v e + DBD modes, experimental limits 
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on half-lives have already been given for 14 out of 34 possible isotopes. The observation 
of 2v e + DBD modes would further constrain the nuclear models employed to study the 
(3(3 decay severely. 

On the other hand, the Oz/ (3(3 decay violates the lepton number conservation and 
is possible in gauge theoretical models beyond the SM as GUTs, Majoron models, R p 
violating SUSY models, lepto quark exchange and compositeness scenario. The aim of all 
the present experimental activities is to observe the Oz/ (3(3 decay. The observation of Oz/ 
c + DBD modes would play a crucial role in discriminating finer issues like the dominance 
of Majorana neutrino mass or right handed currents. The experimental aspects and 
theoretical implications of e + DBD modes have been widely discussed over the past years 
[1-11]. 

The experimental study of (3 (3 decay is usually preferable due to a larger available 
phase space in comparison to e + DBD modes. However, the experimental sensitivity of 
(3~(3~ decay mode gets limited because of the presence of electron background. On the 
other hand, the e + DBD modes are attractive from the experimental point of view due 
to the fact that they can be easily separated from the background contaminations and 
easily detected through coincidence signals from four 7-rays, two 7-rays and one 7-ray for 
(3 + (3 + , (3 + EC and ECEC modes respectively. In the case of the 2z/ ECEC mode, the 

value of 106 Cd is pretty large, 2.782 MeV, but the detection of the + — *0 + transition 
is difficult since only X-rays are emitted. 

In 1955, Winter studied the e + DBD modes of 106 Cd experimentally to explore the 
possibility of distinguishing between the Dirac or Majorana character of the electron 
neutrino [T2] . The 2z/ e + DBD modes were studied theoretically for the first time by Rosen 
and Primakoff pQ. Following the discovery of parity violation in beta decay, there was a 
marked decline in the experimental searches of (3(3 decay in general as both the lepton 
number conservation and the 75 invariance had to be violated for the Oz/ (3(3 decay to occur. 
However, the perception began to change after Vergados showed that e + DBD modes are 
possible as lepton number violating process in gauge theories beyond the SM j2]. Kim 
and Kubodera estimated the half-lives of all the three modes with modified NTMEs and 
non-relativistic phase space factors [T3]. Abad et al. performed similar calculations using 
relativistic Coulomb wave functions 14]. Some other theoretical studies were also done for 
the e + DBD modes [15-18]. The experimental activities on the study of 2z/ e + DBD modes 
were also resumed [19-21]. In the meantime, the QRPA emerged as a successful model 
in explaining the quenching of NTMEs by incorporating the particle-particle part of the 
effective nucleon-nucleon interaction in the proton-neutron channel [22| and the observed 
Ty 2 of several 2z/ (3~ (3~ decay emitters were reproduced successfully [7j. Subsequently, 
the 2z/ e + DBD modes were studied in shell model, QRPA and its extensions, SU^)^ and 
SSDH and pseudo SU(3) [7j. 

Low-background setups using Ge detectors were proposed by Barabash (22] to detect 
the transition of 2z/ ECEC mode to the 0^" excited state. New developments in exper- 
imental setups have led to good limits on the measurement of the 2z/ e + DBD modes of 
nuclei of our interest namely 106 Cd [12,20,24-30] and 108 Cd (21123123 through the direct 
counting experiments. In the mass region A~100, Norman has studied the 2z/ e + DBD 
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modes of 96 Ru [2T] and 102 Pd is also a potential candidate to be studied with Q-value of 
about 1.175 MeV with natural abundance of about 1.02%. With improved sensitivity in 
detection systems of the planned bigger Osaka-OTO experiment [3T] and COBRA [32], it 
is expected that 2v e + DBD modes will be in observable range in the near future. Hence, 
a timely reliable prediction of the half-lives of 96 Ru, 102 Pd, 106 Cd and 108 Cd nuclei will be 
helpful in the ongoing planning of future experimental setups. 

The structure of nuclei in the mass region ^4 pa 100 is quite complex. This mass region 
offers a nice example of shape transitions, i.e. sudden onset of deformation at neutron 
number iV=60. The nuclei are soft vibrators for N<Q0 and quasi rotors for iV> 60. The 
nuclei with neutron number N=60 are transitional nuclei. In this mass region ^4=96-108, 
the smallest and largest quadrupole deformation parameter (3 2 are 0.1580±0.0032 and 
0.2443±0.0030 for 96 Ru and 102 Ru respectively. Further, the pairing of like nucleons plays 
an important role in all (3/3 decay emitters, which are even-Z and even-N nuclei. Thus, 
it is expected that pairing and deformation degrees of freedom will play some crucial 
role in the structure of 96 > 102 Ru, 96 Mo, i02,i06,i08p d and i06,i08 Cd nuclei For the gtudy 

of 2v e + DBD modes of these nuclei, it is desirable to have a framework in which the 
pairing and deformation degrees of freedom are treated on equal footing in its formalism. 
The Projected Hartree-Fock-Bogoliubov (PHFB) model is, in this sense, a sensible choice 
which fulfills these requirements. However in the present version of the PHFB model, it 
is not possible to study the structure of odd-odd nuclei. Hence, the single beta decay 
rates and the distribution of Gamow- Teller strength can not be calculated. On the other 
hand, the study of these processes has implications in the understanding of the role of the 
isoscalar part of the proton-neutron interaction. This is a serious draw back in the present 
formalism of the PHFB model. Notwithstanding, PHFB model has been successfully 
applied to the (3~(3~ decay of many emitters in this mass region, where it was possible to 
describe, in the same context, the lowest excited states of the parent and daughter nuclei, 
as well as their electromangetic transition strengths on one hand, and to reproduce their 
measured double-beta decay rates on the other [33] . 

The aim of nuclear many-body theory is to describe the observed properties of nuclei in 
a coherent framework. The (3(3 decay can be studied in the same framework as many other 
nuclear properties and decays. Experimental studies involving in-beam 7-ray spectroscopy 
concerning the level energies as well as electromagnetic properties has yielded a vast 
amount of data over the past years. Although the availability of data permits a rigorous 
and detailed critique of the ingredients of the microscopic model that seeks to provide a 
description of nuclear (3(3 decay, most of the calculations of 2v e + DBD transition matrix 
elements performed so far but for the work of Barabash et al. [26 a and Suhonen et al. 
[3*1] do not satisfy this criterion. The successful study of 2v e + DBD modes of 106 Cd for 
+ — > + transition together with other observed nuclear properties, like the yrast spectra, 
reduced transition probabilities B(E2: + — > 2 + ), static quadrupole moments Q{2 + ) and 
g -factors g(2 + ) of both parent and daughter nuclei using the PHFB model in conjunction 
with the summation method [35] has motivated us to apply the same framework to study 
the 2v e + DBD modes of 96 Ru, 102 Pd and 108 Cd isotopes. The reason for presenting again 
the results of 106 Cd is that the HFB wave functions are generated with improved accuracy 
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and it is nice to see that the results remain almost unchanged. 

Further, it has been shown that there exists an inverse correlation between the Gamow- 
Teller strength and quadrupole moment jHHl EH]- The PPQQ interaction [38J has two 
terms, associated with the pairing interaction (PP) and the quadrupole-quadrupole (QQ) 
interactions. The former accounts for the sphericity of nucleus, whereas the later increases 
the collectivity in the nuclear intrinsic wave functions and makes the nucleus deformed. 
Hence, the PHFB model in conjunction with the PPQQ interaction is a convenient choice 
to examine the explicit role of deformation on the NTME M 2v . In case of 106 Cd, we have 
already shown that deformation plays an important role in the variation of M 2v vis-a-vis 
changing strength of the QQ part of effective two-body interaction [35J. 

The structure of the present paper is as follows. The theoretical formalism to calculate 
the half-lives of 2v e + DBD modes has been given in a number of reviews [3J [7j and our 
earlier study of 1v e + DBD modes of 106 Cd for the + — > + transition [35J. Hence, we 
briefly outline steps of the above derivations in sect. 2 for clarity in notation. Details 
of the mathematical expressions used to calculate the spectroscopic properties of nuclei 
in the PHFB model have been given by Dixit et al. [39 . In sect. 3, we present results 
to check the reliability of the wave functions of 96 < 102 Ru, 96 Mo, I°2,i06,i08p d and i06,i08 Cd 
nuclei by calculating the mentioned spectroscopic properties and comparing them with the 
available experimental data. The half-lives Ty 2 for the 2z/ e + DBD modes of 96 Ru, 102 Pd, 
106 Cd and 108 Cd nuclei for the + transition are calculated. The role of deformation 
on M 2v is also studied. We present some concluding remarks in sect. 4. 

2 Theoretical framework 

The inverse half-life of the 2u e + DBD mode for the + — > + transition is given by 

T^ 2 (0 + ^0 + )Y 1 = G 2l/ \M 2u \ 2 . (1) 
The G 2l/ is the integrated kinematical factor and the NTME M 2v is expressed as 

M V (0^||ar1|l+)(l+||ar-||0+) 

M 2u - p „ p , {*) 

where 

E = \ (E T - E F ) = l -W . (3) 
The total energy released Wo for different 2z/ e + DBD modes is given by 

W ((3 + (3 + ) = Q p+p+ + 2m e , (4) 
W (f3 + EC) = Q p+EC + e b , (5) 
W (ECEC) = Q E CEC-2m e + e bl + e b2 . (6) 
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The summation over intermediate states is carried out using the summation method 
Pj and the NTME M 2u can be written as 



M 2u = 1- /o+ 



0? , (7) 



where the Gamow- Teller (GT) operator T m has been defined as 

r ™ = ^2 a msT~, (8) 

s 

and 

F m = Y. K -JrD x T mi (9) 

A=0 ^0 

with 

D A r m = [ff, [#, , [#, r m ] ] (A times) . (io) 

When the GT operator commutes with the effective two-body interaction, the Eq. (JJJ) 
can be further simplified to 

m 2v = y. ^B^ijlMI (ii) 

The energy denominator is evaluated as follows. The difference in single particle energies 
of neutrons in the intermediate nucleus and protons in the parent nucleus is mainly due 
to the difference in Coulomb energies. Hence 



Ivi ju) £\P>ir j lit j jn) 



A c - 2E for n„ = n 7r ,l 1/ = l n , j v = j n 

A c - 2E + AEs.o.spitung for n u = n n , l u = l n , j v ^ j w : 

(12) 

where the Coulomb energy difference Ac is given by Bohr and Mottelson 



A 



c 



0.70 

IV3 



(2Z + 1) - 0.76 {{Z + 1) 4/3 - Z 4/3 }] . (13) 



In the case of pseudo SU(3) model [42-44], the energy denominator is a well-defined 
quantity without any free parameter as the GT operator commutes with the two-body 
interaction. The energy denominator has been evaluated exactly for 2v f3~ (3~ HH] and 
2v e + DBD modes [H] in pseudo SU(3) scheme. 

In the present work, we use a Hamiltonian with PPQQ type [3B] of effective two-body 
interaction. Explicitly, the Hamiltonian is written as 

H = H sp + V(P) + Xqg V(QQ), (14) 

where H sp denotes the single particle Hamiltonian. The pairing part of the effective 
two-body interaction V(P) is written as 



V{P) = - SC-l^-^-^otoio^, (15) 

1 a/3 



where a denotes the quantum numbers (nljm). The state a is same as a but with the 
sign of m reversed. The QQ part of the effective interaction V(QQ) is expressed as 

V(QQ) = " (f ) £ £(-ir<a|^|7></%^> 44 ^ (16) 

where 

tf=(iF) ^(0,0). (17) 

The Xgg is an arbitrary adimensional parameter and the final results are obtained by 
setting the x m = 1. The purpose of introducing x m is to study the role of deformation 
by varying the strength of QQ part of the effective two-body interaction. 

The model Hamiltonian given by Eq. is not isospin symmetric. Hence, the 

energy denominator is not as simple as in Eq. (JUJ. However, the violation of isospin 
symmetry for the QQ part of our model Hamiltonian is negligible, as will be evident from 
the parameters of the two-body interaction given later. Further, the violation of isospin 
symmetry for the pairing part of the two-body interaction is presumably small in the 
mass region under study. 

Under these assumptions, the expression to calculate the NTME M 2v of 2u e + DBD 
modes for + — > + transition in the PHFB model is obtained as follows. 

In the PHFB model, states with good angular momentum J are obtained from the 
axially symmetric HFB intrinsic state |$o) with K=0 using the standard projection tech- 
nique jlH] given by 

J D^ (n)R(n)\%)dn, (is) 

where R(Q) and Dq (Q) are the rotation operator and the rotation matrix respectively. 
The axially symmetric HFB intrinsic state |$o) can be written as 

|$o> = Uiuim + v im bt m 4 m )\Q), (19) 

irn 

where the creation operators b\ m and b\ m are defined as 

b%m = Cia,m&\vrn. &Ild 6^ = ^~^ ( 1) 5 Ci aym (X a _ rtl . (20) 
a a 

The results of HFB calculations are summarized by the amplitudes (ui m , Vi m ) and expan- 
sion coefficients Cij >m . 

Finally, one obtains the following expression for NTME Mz v of the 1v e + DBD mode: 



l*o J o) 



(2J+1) 



M> 



2// 



£ 



J f =0, 



00 



o.or t 



I* 



Ji=0 
00 



} 



E + e(n„,l„,j v ) -s(n n ,l 

7V •) J TV ) 



6 



n 



j f =o 
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(21) 



T/ 



where 
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n J = y {det[l + FW(fl)/Wt]}V2 x { de t[l + FM^)/^]} 1 / 2 ^^) sin(0)d0, (22) 
o 

and 



"(z^s-aW*) = {det[l + Fl^W/«U +2 ]} 1/2 X {det[l + F^(0)/^ +2 ]} 1 / 2 . (23) 

The 7r(z/) represents the proton (neutron) of nuclei involved in the 2u e + DBD. The ma- 
trices fz,N and F Z) n{9) are given by 

[/z,Jv]a/3 — J]] CijairriaCijp^mp (l>im a / Ui ma ) 8 ma ^ m ^, (24) 
i 

and 

The calculation of NTME M 2v for the 1v e + DBD mode is carried on as follows. In the 
first step, the matrices [fz,N]a/3 and [F ZtN (9)} a/ s are set up using expressions given by Eqs. 
(24) and (25) respectively. Finally, the required NTME M 2v is calculated in a straight 
forward manner using Eq. (|21|) with 20 gaussian quadrature points in the range (0, 7r). 



3 Results and discussions 

The model space, single particle energies (SPE's) and the effective two-body interaction 
are the same employed in our earlier calculation on 2u e + DBD modes of 106 Cd for the 
+ — > + transition [33]. However, we present a brief discussion of them in the following 
for convenience. The model space consists of lpi/2, 2sx/ 2 , 1^3/2, ld 5 / 2 , 0g7/2, Ogg/2 and 
0/in/2 orbits for protons and neutrons, where we have treated the doubly even nucleus 
76 Sr (N=Z=38) as an inert core. The orbit lpi/2 has been included in the valence space 
to examine the role of the Z=40 proton core vis-a-vis the onset of deformation in highly 
neutron rich isotopes. The set of single particle energies (SPE's) used here are (in MeV) 
e(lp 1/2 )=-0.8, e(0g 9/2 )=0.0, e(ld 5/2 )=5.4, e(25 1/2 )=6.4, e(ld 3 / 2 )=7.9, s(0g 7/2 )=8A and 
e(0h u / 2 )=8.6 for proton and neutrons. This set of SPE's but for the e(0hu/ 2 ), which is 
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slightly lowered, has been employed in a number of successful shell model jlH] as well as 
variational model jU] calculations for nuclear properties in the mass region Af=a 100. 

The strengths of the pairing interaction has been fixed through the relation G p =30/ ^4 
MeV and G n =20/A MeV, which are same as used by Heestand et al. [IB] to explain the 
experimental g{2 + ) data of some even-even Ge, Se, Mo, Ru, Pd, Cd and Te isotopes in 
Greiner's collective model |4^|. The strengths of the like particle components of the QQ 
interaction are taken as: \ vv = Xnn — 0.0105 MeV b~ A , where b is oscillator parameter. 
The strength of proton-neutron (pn) component of the QQ interaction Xpn is varied to fit 
the spectra of 96 ' 102 Ru, 96 Mo, 102 i 106 > 108 p c l and 106 > 108 Cd in agreement with the experimental 
results. To be more specific, we have taken the theoretical spectra to be the optimum 
one if the excitation energy of the 2 + state E 2 + is reproduced as closely as possible to 
the experimental value. Thus for a given model space, SPE's, G p , G n and Xppi we have 
fixed Xpn through the experimentally available energy spectra. We have given the values 
of Xpn in Table 1. These values for the strength of the QQ interaction are comparable to 
those suggested by Arima on the basis of an empirical analysis of the effective two-body 
interactions [201 ■ All the parameters are kept fixed throughout the calculation. 



3.1 The yrast spectra and electromagnetic properties 

In table 1, we have displayed the theoretically calculated and experimentally observed 
values of yrast spectra for J*" =2+, 4+and 6+ states of 96 ' 102 Ru, 96 Mo, 1 02,i06,i08p d and 
106,108 isotopes. The agreement between the experimentally observed (31] and theo- 
retically reproduced E 2 + is quite good. However, it can be noticed that the theoretical 
spectra is more expanded in comparison to the experimental spectra. This can be cor- 
rected to some extent in the PHFB model in conjunction with the VAP prescription [4T] . 
However, our aim is to reproduce properties of the low-lying 2 + state. Hence, we have 
not attempted to invoke the VAP prescription, which will unnecessarily complicate the 
calculations. 

In table 2, we present the calculated as well as the experimentally observed values 
of the reduced B(E2:0 + — > 2 + ) transition probabilities [221, static quadrupole moments 
Q{2 + ) and the gyromagnetic factors g{2 + ) [SH]. In case of B(E2:0 + — > 2 + ), only some 
representative experimental values are tabulated. B(E2:0 + — > 2 + ) results are given for 
effective charges e e ff =0.40, 0.50 and 0.60 in columns 2 to 4, respectively. The experi- 
mentally observed values are displayed in column 5. The calculated values are in excellent 
agreement with the observed B(E2:0 + — > 2 + ) of all the nuclei considered at e e // =0.50 
except for 102 Ru, 102 Pd and 108 Pd, which differ by 0.049, 0.02 and 0.046 e 2 b 2 respectively 
from the experimental lower limits. 

The theoretically calculated Q(2 + ) values are tabulated in columns 6 to 8 of the same 
table 2, along with the experimentally observed Q{2 + ) data in column 9, for the same 
effective charges as used in case of B(E2:0 + — > 2 + ). Again, the agreement between the 
calculated and experimental Q{2 + ) values is quite good in case of 102 Ru, 108 Cd, 106 Pd and 
108 Pd nuclei except for 106 Cd where the difference is 0.2 e b. In case of 96 Ru, 96 Mo and 
102 Pd, although the experimental values have large error bars and a meaningful comparison 
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is difficult, the agreement between calculated and observed values is not satisfactory. The 
g(2 + ) values are calculated with gf =1.0, g"=0.0, g1= g u s =0.60. No experimental 
result for g(2 + ) is available for the isotope 96 Ru. The calculated and experimentally 
observed g{2 + ) values are in excellent agreement for 102 Ru, 102 Pd, 106 Cd and 108 Cd nuclei 
whereas they are off by 0.073, 0.046 and 0.093 nm only for 96 Mo, 106 Pd and 108 Pd isotopes, 
respectively. 

From the above discussions, it is clear that the overall agreement between the cal- 
culated and observed electromagnetic properties is quite good. Hence, the PHFB wave 
functions of 96 ' 102 Ru, 96 Mo, "».io6,i08p d and io6,ios Cd nuclei generate d by fixing Xp n to re- 
produce the yrast spectra are quite reliable. Below, we present the results of NTMEs M 2v 
as well as the half-lives Tffc of 96 Ru, 102 Pd, 106 Cd and 108 Cd for the 0+ -> 0+ transition 
using the same HFB wavef unctions. 

3.2 Results of 2u f3 + f3 + / f3 + EC / ECEC decay 

In table 3, we have compiled the available experimental and theoretical results for 2v 
c+DBD modes of 96 Ru, 102 Pd, 106 Cd and 108 Cd nuclei along with our calculated NTMEs 
M 2v and the corresponding half-lives Ty 2 . The calculated phase space factors were ob- 
tained following the prescription of Doi et al. [I] in the approximation G\ = 1.0, C 2 = 0.0, 
C*3 = 0.0 and Ri t \{e) = R+i(e) + R-i(e) =1.0. The phase space integrals have been eval- 
uated for gp= 1.261 by Doi et al. 0j. However, in heavy nuclei it is more justified to 
use the nuclear matter value of g^ around 1.0. Hence, the theoretical Ty 2 are presented 
both for gA=1.0 and 1.261. 

In case of 96 Ru, the half-life limits Ty 2 of the 2u (3 + EC and 2u ECEC modes for 
the + — > + transition has been investigated by Norman [21] and are of the order of 
10 16 y. The calculated NTMEs M 2v in PHFB and SU(4) CTr model differ by a factor 
of 2 for all the three modes while in QRPA model j^3], the values of NTMEs M 2u are 
larger than the PHFB model values by a factor of 5, approximately. The phase space 
factors for 96 Ru isotope are G 2v ((3 + (3 + ) = 2.516xl0~ 26 y" 1 , G 2u ((3 + EC) = 9.635x 10" 22 
y- 1 and G 2v [ECEC) = 5.385xl0~ 21 y" 1 . The theoretically calculated Tf /2 are of the 
order of 10 26 " 28 y, 10 22 " 23 y and lO 21 " 23 y for 2v P+P+, 2v (3 + EC and 2v ECEC modes 
respectively for g& = 1.261-1.00. 

The e + DBD modes of 102 Pd isotope for + — > + transition has been investigated 
neither experimentally nor theoretically so far. We have used the phase space factors 
G 2v (l3 + EC) = 1.449X10" 30 y" 1 and G 2v {ECEC) = 9.611X lO" 23 y" 1 for 2u f3 + EC and 
2v ECEC modes, respectively. In PHFB model, the predicted Ty 2 of 2v f3 + EC and 2v 
ECEC modes are (2.509-6.344) x 10 32 y and (3.783-9.565) xlO 24 y respectively for g A = 
1.261-1.00. 

We have compiled the available experimental [12,20,24-30] and theoretical results 
[26,34,54-60] for 106 Cd along with our calculated M 2u and corresponding half-life Ty 2 
in table 3. In the case of 106 Cd, the phase factors are G 2u {(3 + (3 + ) = 4.263xl0~ 26 y _1 , 
G 2v (l3 + EC) = 1.570xl0- 21 y^and G 2v {ECEC) = 1.152x 10' 20 y" 1 , respectively. In 
comparison to the theoretically predicted Ty 2 , the present experimental limits for + — >0 + 
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transition of 106 Cd are smaller by a factor of 10 5 7 in case of 2v f3 + (5 + mode but are quite 
close for 2v (3+EC and 2v ECEC modes. The half-life Tfy calculated in PHFB model 
using the summation method differs from all the existing calculations. The presently 
calculated NTME M 2u is smaller than the recently given results in QRPA(WS) model 
of Suhonen and Civitarese [33] by a factor of 2 approximately for all the three modes. 
The theoretical M 2v values of PHFB model and SU(4)o. r [35] again differ by a factor of 
2 approximately for the 2v (3 + EC and 2v ECEC modes. On the other hand, the M 2v 
calculated in our PHFB model is smaller than the values of Hirsch et al. [S3] by a factor 
of 3 approximately in case of 2v j3 + (3 + and 2v ECEC modes while for 2v (3 + EC mode 
the results differ by a factor of 4 approximately. All the rest of the calculations predict 
NTMEs, which are larger than our predicted M 2v approximately by a factor of 7 [571 Iq"%| 
to 10 [2J. The predicted Tfy of 2v /?+/?+, 2v (3+EC and 2v ECEC modes in PHFB 
model are (3.495 - 8.836) xlO 27 y (9.489 - 23.992) xlO 22 y and (1.293 - 3.270) xlO 22 y 
respectively for £4=1.26 and 1.0. 

The 2v ECEC mode of 108 Cd for the + — > + transition has been investigated by 
Georgadze et al. [24] . Kiel et al. [2H] and Danevich et al. |29j . No theoretical calculation 
has been done so far to study the above mentioned mode of 108 Cd isotope. The phase 
space factor of 2v ECEC mode is G 2v (ECEC) = 2.803 xl0~ 26 y" 1 . In PHFB model, the 
calculated half-life T 2 / 2 of the 2v ECEC decay mode is 3.939 xlO 27 y and 9.959 xlO 27 y 
for £4=1-261 and 1.0 respectively. 

The quenching of the nuclear matrix elements seems to be closely related with the 
explicit inclusion of deformation effects, which are absent in the other models. We analyze 
in detail this point below. 

3.3 Deformation effect 

We have investigated the variation of the (Qq), (3 2 and M 2v with respect to the change in 
strength of the QQ interaction x qq to understand the role of deformation on the NTME 
M 2v . Out of several possibilities, we have taken the quadrupole moment of the intrinsic 
state (Qq) (in arbitrary units) and the quadrupole deformation parameter /3 2 as a quanti- 
tative measure of the deformation. The quadrupole moment of the intrinsic states (Qq), 
deformation parameter /3 2 and the NTMEs M 2v for different \ qq are tabulated in table 4. 
The deformation parameter has been calculated with the same effective charge as used in 
the calculation of B(E2:0 + — > 2 + ) transition probabilities. 

It is noticed that the (Qq) as well as f3 2 increases in general as the Xqq is varied from 
to 1.5 except a few anomalies. The intrinsic quadrupole moments show fluctuations in 
case of 96 Ru at Xgg value 0.6. In case of 96 Mo, similar fluctuations are observed at x qq 
equal to 0.4 and 0.6. In case of 102 ' 106 Pd nuclei, the fluctuations occur at x TO =0.4. In all 
cases, it is found that the quadrupole deformation parameter (3 2 follows the same behavior 
as the quadrupole moment of the intrinsic state (Qq) with respect to the change in x qq 
except for 106 Cd isotope. In this case, the (Qq) increases but (3 2 decreases at X(j<j = 0.4. 
Further, there is an anticorrelation between the deformation parameter and the NTME 
M 2v in general but for a few exceptions. 
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To quantify the effect of deformation on M 2v , we define a quantity D 2v as the ratio of 
M 2u at zero deformation (xqq — 0) and full deformation (x qq = 1)- The D 2u is given by 

_ M 2v { Xqq = 0) 

D " - M 2u ( Xqq = iy (26) 

The values of D 2v are 3.13, 3.40, 2.06 and 2.19 for 96 Ru, 102 Pd, 106 Cd and 108 Cd nuclei 
respectively. These values of D 2v suggest that the M 2v is quenched by a factor of 2 to 3.5 
approximately in the mass region 96< A <108 due to deformation effects. 

Given the schematic nature of the PPQQ interaction employed in the present calcula- 
tion, and the fact that many of the nuclei studied are in the transtional region and do not 
display a well defined rotational spectrum, the quenching factors discussed above could 
be considered as a conservative estimate of the uncertainties in the predicted 2v nuclear 
matrix elements. They qualify both the present results and those obtained with other 
models where deformation is not explicitely considered. The uncertainties associated with 
the 0z/ processes would be expected to be far smaller. Future work is expected to clarify 
this point. 



4 Conclusions 

To summarize, we have tested the quality of PHFB wave functions by comparing the 
theoretically calculated results for yrast spectra, reduced B(E2:0 + — > 2 + ) transition 
probabilities, static quadrupole moments Q{2 + ) and g-factors g{2 + ) of 96 ' 102 Ru, 96 Mo, 
I02,i06,i08p ( j an j 106, iosq^ nuc i e j participating in 2v e + DBD modes with the available 
experimental results. The same PHFB wave functions are employed to calculate NTMEs 
M 2v and half-lives Tf /2 of 96 Ru {2v [3+[3+, 2v [3+EC and 2v ECEC modes), 102 Pd (2is 
p+EC and 2v ECEC modes), 106 Cd {2v {3+f3+, 2v (3+EC and 2v ECEC modes) and 
108 Cd [2v ECEC mode) nuclei. It is noticed that the proton-neutron part of the PPQQ 
interaction, which is responsible for triggering deformation in the intrinsic ground state, 
plays an important role in the quenching of M 2v by a factor of approximately 2 to 3.5 in 
the considered mass region 96< A <108. In case of 96 Ru and 106 Cd, we have presented 
and discussed the theoretical results of 2u e + DBD modes in the PHFB model along with 
other available nuclear models for + — > + transition . In case of 102 Pd and 108 Cd, These 
are the first theoretical calculations and in view of growing interests in the study of 2is 
c + DBD modes, these predictions would be helpful in the planning of future experimental 
setups. 
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Table 1. Excitation energies (in MeV) of J 71 " =2+, 4+ and 6+ yrast states of 96 - 102 Ru, 96 Mo, 

102,106,108p d and 106,108 Cd 



Nucleus 


Xpn 




Theo. 


Expt. [51J 


Nucleus 


Xpn 




Theo. 


Expt. [51 


44 Ru 


0.02417 


E 2 + 


0.8323 


0.8326 




0.02557 


E 2 + 


0.7779 


0.7782 






E^+ 


2.1389 


1.51797 








2.0373 


1.6282 








3.8037 


2.1496 






Eq+ 


3.5775 


2.4406 




0.01573 


E 2 + 


0.5551 


0.5565 


44 2 Ru 


0.02054 


E 2 + 


0.4751 


0.4751 






E4+ 


1.6010 


1.2760 






E4+ 


1.4773 


1.1064 






E 6 + 


2.9467 


2.1115 






E G + 


2.8737 


1.8732 




0.01505 


E 2 + 


0.6321 


0.6327 


106 Pr | 
46 ru 


0.01441 


E 2 + 


0.5115 


0.5119 






-E4+ 


1.7298 


1.4939 






-E4+ 


1.4816 


1.2292 






E 6 + 


3.1610 


2.4918 






E G + 


2.7264 


2.0766 




0.01481 


E 2 + 


0.6319 


0.6330 


46 ru 


0.01443 


E 2 + 


0.4336 


0.4339 






E4+ 


1.8072 


1.5084 






E4+ 


1.3126 


1.0482 






E 6 + 


3.3138 


2.5413 






E G + 


2.4826 


1.7712 



Table 2. Comparison of the calculated and experimentally observed reduced transition 
probability B(E2:0 + — ► 2 + ) in e 2 b 2 , static quadrupole moments Q(2 + ) in e b and g- 
factors g(2 + ) in nuclear magneton. Here B(E2) and Q{2 + ) are calculated for effective 
charge e p —l+e e ff and e n = e e ff. The g{2 + ) has been calculated for gf =1.0, g\ =0.0 
and gl = g v s =0.60. 



Nuclei 




^(£2:0+ -> 
Theo. 


2 + ) 
Expt. 52 




Theo. 


Q(2 + ) 


Expt. [23 


Theo. 


S(2 + ) 

Expt. [sni 






e eff 








e eff 












0.40 


0.50 


0.60 




0.40 


0.50 


0.60 








96 Ru 


0.208 


0.261 


0.319 


0.251±0.010 
0.260±0.010 


-0.412 


-0.461 


-0.510 


-0.15±0.27 


0.394 




96 Mo 


0.265 


0.335 


0.413 


0.310±0.047 
0.271±0.005 


-0.466 


-0.524 


-0.582 


-0.20±0.08 
+0.04±0.08 


0.563 


0.419±0.033±0.03^ 


102 Pd 


0.323 


0.410 


0.507 


0.460±0.030 


-0.514 


-0.580 


-0.645 


-0.20±0.20 


0.386 


0.41±0.04 
0.39±0.05 


102 Ru 


0.458 


0.585 


0.726 


0.640±0.006 
0.651±0.016 


-0.613 


-0.692 


-0.771 


-0.57±0.07 
-0.68±0.08 


0.385 


0.371±0.031 


i06 Cd 


0.330 


0.422 


0.525 


0.410±0.020 
0.386±0.05 


-0.518 


-0.586 


-0.654 


-0.28±0.08 


0.372 


0.40±0.10 


106p d 


0.403 


0.515 


0.640 


0.610±0.090 
0.656±0.035 


-0.573 


-0.648 


-0.722 


-0.56±0.08 
-0.51±0.08 


0.465 


0.398±0.021 
0.30±0.06 


IU8 Cd 


0.414 


0.531 


0.661 


0.540±0.011 
0.430±0.020 


-0.581 


-0.657 


-0.734 


-0.45±0.08 


0.361 


0.34±0.09 


108p d 


0.456 


0.584 


0.727 


0.700±0.070 
0.760±0.040 


-0.610 


-0.690 


-0.770 


-0.58±0.04 
-0.51±0.06 


0.483 


0.36±0.03 
0.32±0.03 
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*P.F. Mantica et al. Phy. Rev.C63, 034312 (2001). 
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Table 3. Experimental limits on half-lives Ty 2 (0 + — > + ), theoretically calculated M 2u 
and corresponding T 2 / 2 (0 + -> 0+) for 2v p+/3+, 2u /3+EC and 2v ECEC decay of 96 Ru, 
102 Pd, 106 Cd and 108 Cd nuclei. Half-lives are calculated using g^= (1.261 - 1.0), respec- 
tively. 

Nuclei Decay Experiment Theory 

Mode Ref Tffc (y) Ref. Model \M 2u \ T 2 / 2 (y) 



96 Ru p+p+ 


EH 


> 3.1xl0 16 * 


Present 


PHFB 


0.0537 


(1.378-3.485) xlO 28 








ei 


QRPA 


0.2510 


(6.309-15.950) xlO 26 


(3+EC 


EU 


> 6.7xl0 16 * 


Present 


PHFB 


0.0537 


(3.599-9.100) xlO 23 








E3 


SU(4) (TT 


0.1005 


(1.028-2.598) xlO 23 








El 


QRPA 


0.2694 


(1.430-3.616) xlO 22 


ECEC 




- 


Present 


PHFB 


0.0537 


(0.644-1.628) xlO 23 








m 


SU(4) CTT 


0.1005 


(1.839-4.649) xlO 22 








[54j 


QRPA 


0.2620 


(2.705-6.840) xlO 21 


102 Pd (3+EC 






Present 


PHFB 


0.0524 


(2.509-6.344) xlO 32 


ECEC 






Present 


PHFB 


0.0524 


(3.783-9.565) xlO 24 


10ti Cd /3+/3+ 


|2S1 


>5.0xl0 18 


Present 


PHFB 


0.0819 


(3.495-8.836) xlO 27 




E3 


>2.4xl0 20 ** 


na 


SQRPA(l.b.) 


0.61 


(6.304-15.940) xlO 25 




EB] 


>1.0xl0 19 * 




SQRPA(s.b.) 


0.57 


(7.220-18.260) xlO 25 




E3 


>9.2xl0 17 


El 


QRPA(AWS) 


0.722 


(4.500-11.380) xlO 25 




I2DI 


>2.6xlO m 




QRPA(WS) 


0.166 


(8.513-21.520) xlO 26 




na 


>6xl0 16 


m 


QRPA(WS) 


0.840 


(3.324-8.406) xlO 25 










QRPA(AWS) 


0.780 


(3.856-9.749) xlO 25 








El 


QRPA 


0.218 


(4.936-12.480) xlO 26 










QRPA 




4.940 xlO 25 


/3+EC 


EH 


>1.2xl0 18 


Present 


PHFB 


0.0819 


(9.489-23.992) xlO 22 




E3 


>4.1xl0 20 


m 


SQRPA(l.b.) 


0.61 


(1.712-4.328) xlO 21 




ES 


>0.66xl0 19 * 




SQRPA(s.b.) 


0.57 


(1.960-4.957) xlO 21 




E3 


>2.6xl0 17 


m 


QRPA(AWS) 


0.718 


(1.236-3.124) xlO 21 




EDI 


>5.7xl0 17 * 




QRPA(WS) 


0.168 


(2.257-5.706) xlO 22 








E3 


SU(4) CTT 


0.1947 


(1.680-4.248) xlO 22 








EH1 


RQRPA(AWS) 


0.56 


(2.031-5. 136) xlO 21 










RQRPA(WS) 


0.55 


(2.106-5.324) xlO 21 










QRPA(WS) 


0.84 


(9.027-22.820) xlO 20 










QRPA(AWS) 


0.78 


(1.047-2.647) xlO 21 








El 


QRPA 


0.352 


(5.141-13.000) xlO 21 








EH 


QRPA(WS) 


0.493 


(2.621-6.626) xlO 21 












0.660 


(1.462-3.697) xlO 21 
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Table 3 continued. 





1 'JH 1 

[30J 


>1.UX 1U 


Present 


rnr Jd 


U.Uoly 


{L.2y6-6.2 iu / 


x 1U 




1 OH 1 

[29\ 


>0.oX 1U 




o^KrA^I. D.J 


U.Ol 


(z.ooo-o.oyy / 


x 1U 




1 OO 1 

[28J 


>1.UX 1U 




onn id A /„ u \ 

o^KrA^s.b.J 


U.O ( 


(2.0 < 2-0. fDO y 


nx i n20 
X 1U 




1 O A 1 

123 


>5.8xl0 


123 


l^RPA(AWb) 


0.718 


(1.684-4.258, 


xx i n20 
X 1U 










Llixr J\[ W b ) 


U.loo 


(o.OVo- / . /oU y 


nx 1 f)21 
X 1U 








|59j 


oo-Uri^ i neo j 


U.Zo 


1 1 i n7 o snn^ 


nx i n2i 
X 1U 










OOlJrl^ri/Xp ) 


n i 7 
U. 1 ( 


i''? nn/i 7 c;o.cO 
1^o.uu4- i .oyOy 


X 1U 










SU(4) (7T 


0.1947 


(2.290-5.790) 


xlO 21 










RQRPA(AWS) 


0.56 


(2.768-6.999) 


xlO 20 










RQRPA(WS) 


0.55 


(2.870-7.256) 


xlO 20 








m 


QRPA(WS) 


0.84 


(1.230-3.111) 


xlO 20 










QRPA(AWS) 


0.78 


(1.427-3.608) 


xlO 20 








inn 


QRPA 


0.270 


(1.191-3.011) 


xlO 21 










QRPA(WS) 


0.493 


(3.572-9.031) 


xlO 20 












0.660 


(1.993-5.039) 


xlO 20 


108 Cd ECEC 


1291 


>4.1xl0 17 


Present 


PHFB 


0.0952 


(3.939-9.959) 


xlO 27 



[28J >1.0xl0 18 

[24J >4.1xl0 17 

* and ** denote half-life limit for 0u + 2v and Oz/ + 1v + Oz/M modes respectively. 
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Table 4.Effect of the variation in Xgg on {Qp, fa and NTMEs M 2v . 





Xqq 


0.00 


0.20 


0.40 


0.60 


0.80 


0.90 


0.95 


1.00 


1.05 


1.10 


1.20 


1.30 


1 


96 Ru 


{QD 


0.0 


0.006 


0.214 


0.144 


23.854 


30.351 


32.485 


34.473 


36.42 


38.15 


66.54 


70.05 


7; 




09 


0.0 


0.046 


0.097 


0.098 


0.112 


0.140 


0.151 


0.161 


0.171 


0.180 


0.295 


0.317 


o 


96 Mo 


{QD 


0.0 


0.695 


0.211 


0.477 


22.464 


31.82 


37.02 


41.73 


45.15 


48.15 


61.43 


65.44 


6i 




fa 


0.0 


0.091 


0.093 


0.093 


0.106 


0.149 


0.174 


0.191 


0.210 


0.224 


0.268 


0.281 







M 2u 


0.168 


0.154 


0.152 


0.153 


0.093 


0.072 


0.067 


0.054 


0.036 


0.024 


0.076 


0.049 





102p d 


{QD 


0.0 


0.252 


0.081 


1.080 


1.839 


36.08 


42.32 


45.47 


47.91 


49.63 


52.67 


56.37 


8- 




09 


0.0 


0.085 


0.046 


0.090 


0.092 


0.149 


0.172 


0.185 


0.196 


0.203 


0.217 


0.234 





102 Ru 


{QD 


0.0 


0.028 


0.123 


0.492 


38.04 


49.81 


53.57 


56.51 


58.87 


61.19 


66.93 


88.17 


8; 




fa 


0.0 


0.036 


0.068 


0.092 


0.159 


0.204 


0.220 


0.232 


0.242 


0.252 


0.279 


0.362 









0.178 


0.203 


0.208 


0.215 


0.135 


0.092 


0.072 


0.052 


0.039 


0.027 


0.021 


0.0001 





106 Cd 


{QD 


0.0 


0.008 


0.031 


0.128 


0.510 


32.65 


40.83 


47.14 


55.89 


62.54 


73.53 


83.65 


91 




fa 


0.0 


0.007 


0.003 


0.035 


0.073 


0.127 


0.152 


0.176 


0.211 


0.243 


0.299 


0.325 





106p d 


{QD 


0.0 


0.022 


0.079 


0.192 


0.897 


39.54 


46.88 


52.12 


56.13 


59.31 


66.23 


73.87 


7! 




fa 


0.0 


0.016 


0.042 


0.064 


0.094 


0.158 


0.183 


0.203 


0.216 


0.227 


0.254 


0.290 







M 2v 


0.169 


0.164 


0.162 


0.166 


0.170 


0.127 


0.095 


0.082 


0.084 


0.066 


0.041 


0.027 





108 Cd 


{QD 


0.0 


0.015 


0.047 


0.132 


0.488 


35.73 


43.87 


53.29 


67.98 


76.17 


79.26 


80.61 


8< 




fa 


0.0 


0.001 


0.013 


0.039 


0.076 


0.134 


0.161 


0.195 


0.258 


0.299 


0.311 


0.316 





108p d 


{QD 


0.0 


0.047 


0.100 


0.267 


28.21 


44.76 


50.99 


55.88 


59.63 


62.85 


70.17 


76.61 


8; 




fa 


0.0 


0.046 


0.051 


0.076 


0.125 


0.174 


0.196 


0.213 


0.225 


0.236 


0.267 


0.299 







M 2v 


0.208 


0.199 


0.203 


0.204 


0.183 


0.120 


0.093 


0.095 


0.084 


0.052 


0.029 


0.020 
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